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In recent years a unified phenomenological picture for tlie hole doped high- Tc 
cuprates has emerged for a spin and charge spectroscopy. Spectral anomalies have 
been interpreted as evidence of charge carrier coupling to a collective spin excitation 
present in the optical conductivity, in ARPES (angular resolved photoemission), 
and in tunneling data. These anomalies can be used to derive an approximate 
picture of a charge carrier-exchange boson interaction spectral density I'^xi^) which 
is then be used within an extended Eliashberg formalism to analyze normal and 
superconducting properties of optimally doped and overdoped cuprates. This paper 
reviews recent developments and demonstrates the sometimes astonishing agreement 
between experiment and theoretical prediction. 
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INTRODUCTION 

The standard Eliashberg equations 
were derived for superconductors with an 
energy gap of s-wave symmetry and the 
electron-phonon interaction as the pairing in- 
teraction. This type of interaction allows 
the application of Migdal's theorem which 
states that vertex corrections in the electron- 
phonon interaction can be neglected to or- 
der uj£,/ep, with lud the Debye energy and 
6f the Fermi energy. On the other hand, 
it is now widely accepted that the high Tc 
cuprates have an energy gap of d^i-yi sym- 
metry [§, H, H, ^ 0, H and there is still no 
consensus as to the microscopic mechanism 
leading to Cooper pairs in these materials. 

In Eliashberg theory a given superconduc- 
tor is characterized by the Eliashberg func- 
tion a'^F{uj) which describes the exchange 
of a phonon by two electrons at the Fermi 
surface and by the Coulomb potential /i*. 
These are the kernels in the two non-linear 



coupled Ehashberg equations. One equation, 
which is referred to as the renormalization 
channel, describes the effect of the electron- 
phonon interaction on normal-state proper- 
ties modified further by the onset of super- 
conductivity. The second equation, referred 
to as the pairing channel, deals with the en- 
ergy gap directly and is identically zero in 
the normal-state. When reliable tunneling 
data is available for the quasiparticle density- 
of-states Nqp{uj), for instance, the procedure 
can be inverted |]^, [T^ to get from Nqp 
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the kernels a^Fiuj) and /i*. In principle, it 
should also be possible to get the same in- 
formation from infrared data |Tl|, |12|, ^ al- 
though in conventional systems this has not 
been widely done while tunneling has. Once 
the kernels q^F{(jj) and /x* are known, the fi- 
nite temperature Eliashberg equations can be 
solved numerically to obtain superconducting 
properties. 

In principle, the Eliashberg equations can 
easily be generalized to include (i-wave sym- 
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metry of the energy gap by an appropriate 
extension (to include a dependence on orien- 
tation of the electron momenta) of the charge 
carrier-exchange boson interaction spectral 
density {a'^F{uj) in case of the electron- 
phonon interaction) which contains all the 
relevant information about the coupling of 
the charge carriers to the exchange bosons. 
As the microscopic mechanism leading to su- 
perconductivity is not yet known, informa- 
tion on the charge carrier-exchange boson in- 
teraction spectral density (denoted Px{^) 
throughout this paper) is to be obtained by a 
fit to appropriate data sets using phenomeno- 
logical models. Such a procedure can yield a 
first approximation to a complete description 
even in cases when an equivalent to Migdal's 
theorem is not applicable and vertex correc- 
tions are not entirely negligible. 

Unfortunately, the well established inver- 
sion techniques which allowed one to de- 
termine a'^F{u!) from tunneling experiments 
p!0| have, so far, not been extended to the 
cuprates and, therefore, phenomenological 
models had to be developed for Pxi^)- Oiie 
such phenomenological model has been intro 
duced by Schachinger et al 



in optimally doped YBa2Cu306.95 (YBCO) 
pO| , 21], the similar peak observed in the 
electronic thermal conductivity |2^, and the 
temperature dependence of the penetration 
depth in nominally pure YBCO samples and 
in YBCO samples with Zn or Ni impurities 
I^ . Nevertheless, similar results would have 
been achieved using the Marginal Fermi Liq- 
uid (MFL) model assuming a d-wave gap 
together with a low frequency cutoff to de- 
scribe the charge carrier-exchange boson in- 
teraction spectral density. 



A remarkable step forward in the devel- 
opment of phenomenological models was pro- 
vided by the work of Marsiglio et al. 



21 who 



were able to show analytically that there ex- 
ists a simple, approximate formula which re- 
lates a^F[uj) to the normal state optical con- 
ductivity (j{uj) via the second derivative of the 
real part of uja~^{uj). This established the ba- 
sis for a spectroscopy which allows the mea- 
surement of the spectral density a^Fiu) di- 
rectly from optical data. This result was then 
extended to the superconducting state of d- 



1^ O n!6l and wave superconductors by Carbotte et al. |^ 



was reviewed by Schachinger and Schiirrer 
r7| . This model is a purely electronic model 



and describes the feedback effect the super- 
conducting state has on /^x(a;). The au- 
thors used, for definiteness, the spin fluctua- 
tion model introduced by Pines and cowork- 
ers |19| in their Nearly Antiferromag- 
netic Fermi Liquid (NAFFL) model. The 
feedback effect caused by superconductiv- 
ity is described by introducing a low en- 
ergy gap in Pxiy^) (low frequency cutoff) 
which opens up as the temperature is lowered 
through the critical temperature T^. This 
gap shows the same temperature dependence 
and size as the superconducting energy gap. 
Within this model it was possible to de- 
scribe consistently the temperature depen- 
dence of the microwave conductivity with 
its pronounced peak around 40 K observed 



who explore the relationship between spec- 
tral density and W{uj) at low temperatures 
in the superconducting state. They conclude 
that the relationship is not at all as direct, 
but, even though more complicated, it re- 
mains simple enough to be very useful al- 
though more approximate. (A similar pro- 
cedure was also suggested by Munzar et al. 



It will be the purpose of this paper to re- 
view the application of this technique to var- 
ious cuprates in some detail. Thus, the pa- 
per establishes in section two the formalism, 
section three discusses its application to opti- 
mally doped YBCO. Other cuprates are also 
investigated within the same context, and fi- 
nally, in section four a summary is presented. 
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FORMALISM 

The Normal State Optical Conductivity 

The optical conductivity is related to the 
current-current correlation function. The 
paramagnetic part of the response function 
on the imaginary frequency axis is given by 

0, ^, ^ §n|, g 



tr < G(k, iuj„ 



(1) 



where G(k, iuo^) is a matrix Green's func- 
tion in the Nambu formalism p2|, iujm = 
i7i'T{2m+l),m = 0, ±1, ±2, ... is the fermion 
and ii'n = 2in7iT, n = 0, ±1, ±2, ... is the bo- 
son Matsubara frequency; T is the tempera- 
ture and Vx the component of the electron 
velocity in x-direction. The factors preced- 
ing the summations include the total number 
of atoms in the crystal, A^, and the inverse 
temperature, /3 = l/fc^T. 

The optical conductivity is related to the 
response function through 



a{uj) = -n(cj + iO+). 

UJ 



(2) 



After analytical continuation to the real fre- 
quency axis and using the usual procedure 



N ^ 



de N{e) 



we arrive at a general expression for the op- 
tical conductivity a'{uj): 



a'iu;) 



du tanh 



U + UJ 

2T 



xS-\T,uj,u) 



tanh 







— tanh I — 



V2T7J 



U + UJ 

2T 

S-\T,uj,u)\,{3) 



which has been given by Lee et al. p8[ (with 
the factor ne^ /m suppressed). Here, N{e) = 
Niep) = iV(0), N{0)e^vl = QI/Ati = ne^/m; 
Qp is the plasma frequency, e the charge on 
the electron, m its mass, and n the electron 
density per unit volume. In Eq. (H) 



-S(T,z/) 



(4) 



with the self energy S(r, cj) related to the 
electron-phonon spectral density by 

S(T,^) = - jdza'Fiz) + 



, , 1 .UJ - z 

-lb { - + t 

^ ' 2 27iT 



(5) 



where ip is the digamma function. In Eq. (1) 
nt~^ = l/{2Timp) gives the impurity contribu- 
tion to the electronic scattering, (rj^p is the 
impurity scattering time.) 

At zero temperature these expressions for 



the conductivity reduce to a simple form p4, 
3|: 



a(uj) = — — / du 



1 



An UJ J UJ — S(z/) — S(a; — u) 





and the self energy 



(6) 



oo 

S(a;) = jdna^F{n)\n 



-in jdVta^F{Vt) 



for the electron-phonon interaction. In this 
form Marsiglio et al. [0] were able to show 



analytically that a remarkably simple formula 
could be used to establish an approximate 
but very useful relationship between (7{uj) and 
a'^F{uj). The observation was also backed 
up by detailed numerical work. We begin by 
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FIG. 1: The solid curve is a^F^uj) vs. to for 
Pb. The other curves are the function W{uj) 
according to Eq. (^) obtained from the normal- 
state conductivity optical scattering rate t~p^(lli) 
at various temperatures. See Marsiglio et al. 



defining an optical scattering rate t^{uj) as 



24, 25 



-^3?e- 



Top{oj) 4:71 a (to 



3fJe- 



(7) 



which is routinely obtained in optical experi- 
ments. We then define an auxiliary function 



24, 25 



W{uj) 



1 rf2 



27r dto'^ 1 Top{uj) 



UJ 



Marsiglio et al. have shown that in cer- 
tain circumstances 



from Eq. (H) and computing W^u) defined 
by Eq. (H). Such calculations were performed 
at two temperatures, namely T = 1 K (dot- 
ted) and T = UK (dash-dotted). Within 
the energy range corresponding to the range 
of a'^F^u) the dotted curve for W{uj) is re- 
markably close to the solid curve for a'^F{uj) 
and therefore W{u) gives an accurate mea- 
surement of the absolute value as well as the 
frequency dependence of the spectral density. 
As the temperature is increased this is no 
longer the case although some rough corre- 
spondence remains which provides a quali- 
tative similarity between the two quantities 
which could still be exploited to get a rough 
first measure of the spectral density in cases 
where low temperature data are not available. 
We point out that even at T = IK there 
are negative tails in W{uj) above the maxi- 
mum phonon cutoff which are not in a'^F{u). 
This is expected since W{u) and a^F(a;) are 
not the same quantities. In fact, it is indeed 
remarkable that they should correspond so 
closely below the phonon cutoff energy. This 
close correspondence can be exploited to get 
a good first measure of the spectral density 
a'^F{uj) from infrared data. In principle, one 
should use a first iteration for a^F^uj) ob- 
tained from the second derivative of the con- 
ductivity defining W{lj) [Eq. (H)], to calculate 
from it <j{uj) based on Eqs. @) to (|^) and keep 
iterating until an exact correspondence be- 
tween calculated and measured cr{uj) results 
has been achieved. In any of the applications 
so far this has not been attempted because of 
the many uncertainties that remain. 



(9) 



which serves as a basis for a spectroscopy 
which allows the measurement of the spectral 
density a'^F^u) directly from optical data. In 
Fig. ID we show theoretical results for W{uj) 
at two temperatures based on the case of 
Pb. The solid curve is the Pb a'^F{u!) ob- 
tained from tunneling data. The other two 
curves were obtained by calculating cr^u) 



The Superconducting State 

In as much as BCS theory applies, any ef- 
fective interaction between two electrons at 
the Fermi surface, which is attractive, will 
lead to superconductivity. This can arise 
from the electron-phonon interaction through 
the polarization of the system of ions. In this 
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case we can describe the polarization process 
as due to the exchange of a phonon between 
a pair of charge carriers. An obvious ex- 
tension is to ask: could we exchange some 
other excitation? In the Nearly Antiferro- 
magnetic Fermi Liquid (NAFFL) model of 
Pines and coworkers |T8|, it is envisaged 
that spin fluctuations replace the phonons. 
The basic formalism for dealing with this new 
situation are the Eliashberg equations but 
now the k, k' anisotropy in momentum of 
the kernel a1^^,F{uj) needs to be taken into 
account so that the resulting superconduct- 
ing state exhibits d-wave symmetry to ac- 
cord with the experimental observation. In 
the renormalization channel we take for sim- 
plicity only the isotropic contribution from 
the electron-spin fluctuation exchange writ- 
ten as /^x(a;) where P is to denote a spin- 
charge exchange coupling constant and x{^) 
is the spin susceptibility. To get a rf-wave 
gap we use a separable interaction of the form 
cos{2e)gPx{uj) cos(2^') with 6 and 6' the di- 
rection of the initial (k) and final (k') mo- 
mentum which, for simplicity, we pin on the 
Fermi surface although in the NAFFL the en- 
tire Brillouin zone is averaged over i.e.: is 
not pinned on the Fermi surface. While, for 
simplicity, we have assumed the same form 
Px{^) to hold in the pairing as in the renor- 
malization channel we have introduced a nu- 
merical factor g to account for the fact that 
the projection of the general spectral density 
will in general be different in the two chan- 
nels. The repulsive effective Coulomb inter- 
action fi^ j^, is isotropic in an isotropic s-wave 
formalism. The same holds for the 'Hub- 
bard' U which is also assumed to be large 
and isotropic. Thus, the effective Coulomb 
potential is not expected to have a numeri- 
cally large ci-wave symmetric part and there- 
fore does not contribute to the pairing chan- 
nel in a generalized c?-wave formulation of the 
Eliashberg equations. 

Within these simplifying assumptions the 
Eliashberg equations need first to be written 



on the imaginary Matsubara frequency axis. 
They take on the following form |3^: 



A{iun;9) = ^7rT^cos(2^)A( 



m — n] 




'10a) 



for the renormalized pairing potential 
A{iLJn] 0), and 

u{iujn) = uin + ttT A(m — n) 




(10b) 

for the renormalized frequencies ui{iuJn)- 
Here, (■ ■ ■ ) denotes the angular average over 
6. The quantity A(m — n) has the usual form 



A(m — n) 



UJr, 



:■ (10c) 



As written, Eqs. (0) do not depend on im- 
purity scattering. To include this possibility, 
we need to add into the right hand side of Eq. 
( |10bD a term of the form 

r+^^ — „„,.\"^_„,. , (11) 



TT 



where F"*" is proportional to the impurity 
concentration and c is related to the elec- 
tron phase shift for scattering off the impu- 
rity. For unitary scattering, c is equal to zero 
while c ^ oo gives the Born approximation, 
i.e.: the weak scattering limit. In this case 
the entire impurity term reduces to the form 
'nt^^[iiOn) with c absorbed into t"*". To com- 
plete the specification of Eq. (pT]) , we have 



AiiuJr. 



(12a) 
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and rectly from the numerical solution on the 

imaginary frequency axis, i.e.: from A(ic(j„; 6) 




and oj{iUn), real frequency axis solutions are 
+ A'^iiuj -6)1 needed to calculate the optical conductivity. 

These equations for A(z/ + i5; 9) and u){v + i5) 
with 5 infinitesimal are more complicated and 



l^t^..., 

(12b) 

While certain quantities, such as the pen- ^^^^^^^ ^^^^ 0,00: 

etration depth, can be obtained quite di- 



A(z/ + i5] 9) = nTg cos(26') [A(z/ — ium) + A(z/ + ium)] 



m=0 

K{iuJrn;9')cos{29' 

X 



\J UJ^{iuJm) + /\?{iuJm] 9') 

oo 

+iTig j dz cos(26')/\(z) [n{z) + f{z - v)] x 

-oo 

K{v - z + i5;9')cos{29') 



'u'^iiy -z + t6)- A2(z/ -z + t6; 9') 

for the pairing channel and 

oo 

u){v + i5) = u + inT [A(z/ — icOm) — A(z/ + icom)] 



m=0 

/ 



+ 



:i3a) 



'uj'^{iujm) + A2(iu;„; 9' 

— oo 

I 

for the renormalization channel. Thermal Bose and Fermi distribution n{z) and f{z), 
factors appear in these equations through the respectively. Furthermore, the abbreviations: 
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troduced in the usual way as u^u + iS) = 
vZiy)^ one finds for tlie superconducting gap 



Div) 



(14b) 



(14c) 

liave been used. Eqs. ([13|) are a set of non- 
linear coupled equations for the renormalized 
pairing potential A(z/ + i6; 9) and the renor- 
malized frequencies ujiv + id) with the super- 
conducting gap given by 

A(. + ,M) = .toM), (15) 

UJ^U + lo) 

or, if the renormalization function Z{v) is in- 



Aiv + i5: 



A(i^ + z5;i 



(16) 



These equations are a minimum set and go 
beyond a BCS approach and include the in- 
elastic scattering known to be strong in the 
cuprate superconductors. 

From solutions of the generalized Eliash- 
berg equations we can construct the Green's 
function in Eq. ([l|) analytically continued to 
the real frequency axis VL. In this formulation 
the expression for the in-plane conductivity 
(Jab [T, fi) involves further averaging over an- 
gles which needs to be done numerically. We 
find the following result after further manip- 
ulations and rearrangements: 



i e'^N{0)vl 
n 2 



, , / u \l-N{i^;9)N{u + n;e)-Pii^;9)P{i^ + n;^ 

X ( / dz/tanh — - — — — — — — 

\2TJ E{v-9) + E{v + vt-e) 

'0 



+ /d.tanh(^) 



v + 1 - Ar*(z/; 9)N%v + VL]9)- P*{u; 9)P*{u + Q; i 



E*(u:9)+E*(u + n: 



X 



tanh ( — I — tanh ( —=] 
2T / \2TJ 



1 + N*{v] 9)N{v + n]9) + P*{u] 9)P{u + fi; 



E{v + n-9) - E*{v- 9) 

v + r 1 - N*{v] 9)N*{iy + n;9)- P*{iy; 9)P*{v + fi; 9) 



EHv- 9) + E*(u + Q: 



1 + N^u; 9)N{v + Vt-9) + P^(z/; 9)P{v + fi; 
E(iy + n-9) - E*(iy- 9) 



(17a) 
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with 



and 



N{uj-e) 



E{iu;e) 

E(uj-e) 



(17b) 

(17c) 
(17d) 



In the above, the star refers to the complex 
conjugate. This set of equations is vahd for 
the real and imaginary part of the conduc- 
tivity function of frequency fl. It con- 
tains only the paramagnetic contribution to 
the conductivity but this is fine since we have 
found that the diamagnetic contribution is 
small in the case considered here. 



The out-of-plane conductivity <7c(T, fl) at 
temperature T and frequency Q is related 
to the current- current correlation function 
UciTjiUn) at the boson Matsubara frequency 
iun analytically continued to real frequency 
Q, and to the c-axis kinetic energy {He) 



via 



[ne(T, 



with d the distance between planes in c- 
direction. In terms of the in-plane thermo- 
dynamic Green's function G{k, iujn) and for 
coherent hopping t_L(k) perpendicular to the 
Cu02 planes 



ne(T,zz/„) = 2(erf)2T^^ti(k)tr|foG(k,icj„)foG(k,icj„ + iz/„ 



•^m k 



and 



{H,) = 2TJ2Y1 ^i(^) {^^^(k, 2^„)f3G'(k, tiu, 

I 



(19a) 
(19b) 



l^m k 



In Eqs. ( p!9D the 2x2 Nambu Green's function 
G(k, ium) describes the in-plane dynamics of 
the charge carriers with momentum k in the 
two dimensional Cu02 plane Brillouin zone 
and is given by 



G{k,iuJn) 



iu{iUn)fo + CkT-3 + Ak(iu;„)'ri 



—Uj'^iiuJr, 



(20) 

where the f's are Pauli 2x2 matrices, 
is the band energy of the charge carriers as a 
function of their momentum k, A]^{iujn) is the 
renormalized pairing potential and iuiiun) 
the renormalized Matsubara frequency. In 
our model these quantities are determined as 
solutions of Eliashberg equations (p!0|). 

In Eqs. (^) the out-of-plane matrix ele- 



ment t_L(k) can depend on the in-plane mo- 
mentum k. Models have been summarized 
recently by Sandeman and Schofield who 



refer to previous literature |3^, A pos- 
sible choice is t±(k) = a constant. But, 
consideration of the chemistry of the Cu02 
plane and of the overlap of one plane with 
the next, suggests a form t^{\s) = cos^(2^) 
where 6 is the angle of k in the two dimen- 
sional Cu02 Brillouin zone for the plane mo- 
tion. This matrix element eliminates entirely 
contributions from nodal quasiparticles to the 
c-axis motion. 

For incoherent impurity induced c-axis 
charge transfer Eqs. (|19D are to be modified. 
After an impurity configuration average one 
obtains 
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n,(T,^z/„) = 2(erf)2r^^Vk2k,tr|foG(k,zcu^)foG(k',^cu„ + ^z/„)| 



m k,k' 



m = 2TJ] J]l^4,tr{f3G(k,za;^)f3G(k', 



(21a) 
(21b) 



m k,k' 



with V^]^, the average of the square of the 
impurity potential If the impurity potential 
was taken to conserve momentum, which it 
does not, we would recover Eqs. (|T9D. Var- 
ious models could be taken for V^-^^,. Here 
we use a form introduced by Kim and 
Hirschfeld et al. IHll 



*^k,k' 



+ |14pcos(2^)cos(20'), (22) 



with 9 and 9' the directions of k and k' re- 
spectively. 

After analytic continuation to the real Vt- 
axis the real part of the incoherent conductiv- 
ity along the c-axis is given by (normalized to 
its normal state value aicn) P] : 



(Jic(fi) 



+ 



(23) 



APPLICATION TO HIGH-Tc 
CUPRATES 

The Compound YBa2Cu306+5 

The Normal State Infrared Conductivity 

We begin with a discussion of the normal- 
state scattering of the charge carriers off spin 
fluctuations. For the spin fluctuation spectral 
density Pxi^) ^6 take a very simple model 



(MMP). We define a single characteristic spin 
fluctuation frequency usf and take 



Uj/USF 



(24) 



SF 



where P is a coupling constant which can 
be fit to normal-state infrared data based 
on Eq. (0) with J^x(u;) playing the role of 
a^F{uj) in Eq. (^ for the self energy. 

In Fig. H we show our result for t^{ijj) 
related to the conductivity by Eq. (|^) for 
an optimally doped, twinned YBa2Cu306.95 
(YBCO) single crystal with = 92.4 K at 
temperature T = 95 K. The solid curve is 
the data of Basov et al. WM- The dashed 



motivated in the work of Millis et al. [ITS 



and dotted curves are our best fits for ujsf = 
10 meV and 30 meV, respectively, with P ad- 
justed to get the correct absolute value of the 
scattering rate at T = 95 K and low ener- 
gies UJ. We see that both values of ujsf do 
not give equally satisfactory fits to the data. 
The dash-dotted curve, however, fits the data 
well and corresponds to usf = 20meV. This 
fit provides us with a model Pxif^) valid 
for the normal-state of YBCO. A plot of 
this function is shown in Fig. ^ as the gray 
solid squares. Also shown in this figure are 
two sets of theoretical results based on the 
Pxi'^) with ujsF = 20meV which serve to 
illustrate the inversion technique. The ex- 
perimental data on r^ioj) gives the model 
Pxiyj) spectrum. Next this model spec- 
trum is used in the normal-state conductiv- 
ity Eq. d^) and two temperatures are consid- 
ered, namely T = 95 K (dotted curve) and 
T = lOK (solid curve). Except for some 
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FIG. 2: The normal-state optical scattering rate 
ropi(Lj) vs uj for YBCO with a = 92.4 K ob- 
tained from the work of Basov et al. [42| (solid 
line) at a temperature of 95 K. The dash-dotted 
curve from theory based on Eq. (^) with an MMP 
model spectral density using a spin fluctuation 
frequency ujsf = 20meV gives good agreement 
while the other choices of 30meV (dotted line) 
or lOmeV (dashed line) do not. 



numerical noise it is clear that when W{uj) 
(solid curve) is computed from Eq. (|^) which 
involves a second derivative of our calculated 
data for t~^{uj) computed from the conduc- 
tivity according to Eq. (^, that the low tem- 
perature data (T = 10 K) gives a remarkable 
accurate picture of the input Pxif^) function 
and that to a very good approximation Wiu) 
is the same as Px{'^) this case. If the tem- 
perature is increased to 95 K the match be- 
tween W{(jj) (dotted curve) and Pxi}^) is not 
quite as good. This shows that our inversion 
technique summarized in Eqs. (|^ and (||) is 
best when applied at low temperatures. The 
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FIG. 3: Comparison of the spectral density 
l'^x{u}) (gray solid squares) in the MMP model 
jlSll {ojsF = 20meV with the function W{uj) de- 
fined in Eq. (|8|) from the normal-state conduc- 
tivity scattering rate Eq. (^. The solid curve is 
at T = lOK while the dotted curve is at 95 K. 

resultant inverted Pxi}^) gets smeared some- 
what if high temperatures are used instead. 

The Superconducting State, I'^xi^) o^^^d the 
Infrared Conductivity 

Going to the superconducting state re- 
quires the solution of the Eliashberg equa- 
tions ([T3| ) and evaluation of formula (|1^) for 
the optical conductivity. This is much more 
complicated than the corresponding normal- 
state analysis. Also, it is critical to under- 
stand that since we are dealing with a highly 
correlated system and the excitations we are 
exchanging in our spectral density are within 
the electronic system itself, they could be pro- 
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Energy (meV) 

FIG. 4: Experimental results of the mag- 
netic susceptibility at Q = (vr, vr) in a sam- 
ple of YBa2Cu306.92- The imaginary part of 
X is measured as a function of energy. Open 
squares are results at T = 100 K and solid circles 
are at T = 5K (in the superconducting state, 
Tc = 91 K for this sample). The energy of the 
spin resonance Ej. = 41 meV. Adapted from . 

foundly modified by the onset of the tran- 
sition. In Fig. H we reproduce spin polar- 
ized inelastic neutron scattering results of 
the spin fluctuations measured in a sample 
of YBa2Cu3 06.92 by Bourges et al. |Q at 
two temperatures namely T = 100 K (open 
squares) and T = 5 K (solid circles) which 
show the formation of the 41 meV spin reso- 
nance in the superconducting state. At 100 K 
the measured spectrum looks very much like 
the simple spectrum used in our analysis of 
the optical data in YBCO (Fig. ||) but this 
simple form is profoundly modified in the su- 
perconducting state with the imaginary part 
of the magnetic susceptibility depressed at 
small uj and the formation of a sharp peak 
around = 41 meV. The possibility of a 
change in Px{^) brought about by the on- 
set of superconductivity must be included in 
our analysis of the optical data in the super- 
conducting state. A question we can immedi- 



ately ask is: is the 41 meV peak seen in opti- 
mally doped YBCO in neutron scattering also 
seen in the superconducting state optical con- 
ductivity? To perform the necessary analysis 
several modifications of what has been done 
so far need to be considered. The observa- 
tion that, in the normal-state the well defined 
function W{u!) given in Eq. (^ which is eas- 
ily accessible when the conductivity cr(u;) is 
known, can be identified to a good approxi- 
mation with J^x(u;) may not hold in the su- 
perconducting state of a (i-wave superconduc- 
tor. 

It turns out from extensive calculations of 
the superconducting state conductivity based 
on Eqs. {^i7\) by Schachinger and Carbotte 
|]4^ , ^] that a simple modification of the rule 

can be found which applies approximately in 
the superconducting state in the resonance 
region, namely 

with UJ2 shifted by Aq(T), the energy gap 
at temperature T, when compared with uji. 
This rule, while not exact, is nevertheless suf- 
ficiently accurate to make it useful in obtain- 
ing first qualitative information on the mag- 
nitude and frequency dependence of the un- 
derlying charge carrier-exchange boson inter- 
action spectral density from optical data. We 
expect the changes in going to the supercon- 
ducting state to be the growth of the reso- 
nance, while at higher energies there should 
be no change in Px{^) from its normal state 
value. A detailed comparison of Px{^) 
W{uj) in the superconducting state will be 
given later when we relate additional struc- 
tures in W{u!) not in Px{^) to structures in 
the superconducting quasiparticle density of 
states which introduce distortions in W^u) 
as compared to the underlying spectral den- 
sity. In Fig. ^ we show results for the case 
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FIG. 5: The model charge carrier spin ex- 
citation spectral density I'^xi^) ^oi T = lOK 
(gray solid squares) constructed from conductiv- 
ity data for optimally doped YBCO. The dashed 
line which follows the gray solid squares faith- 
fully, except for negative oscillations just beyond 
the spin resonance around 68meV, is W{uj)/2 
obtained from our model I'^xi^) (displaced by 
the gap energy Aq = 27meV in the figure). The 
solid line is the coupling to the resonance found 
directly from experiment. 



of YBCO ^ H at T = 10 K. The solid 
line was obtained directly from differentia- 
tion of experimental data on the optical scat- 
tering rate using the superconducting state 
conductivity cr^uj) in Eq. (^ and the defini- 
tion (^ of W{uj). It is limited to the reso- 
nance region. It shows that W{uj)/2 clearly 
has a peak corresponding to the peak seen 
in the spin susceptibility of Fig. § measured 
in neutron scattering. Of course, Px{^) i^^" 
eludes the coupling constant P between elec- 
trons and spin fluctuations and is not strictly 
the imaginary part of the spin susceptibil- 
ity. These two functions are not the same 
but we do know that an optical resonance 
peak does appear when the neutron peak is 
observed. They fall at the same frequency 
Er and look similar in other aspects. It is 
important to emphasize that this solid curve 
comes directly from the data on (t{uj) and can 
be interpreted as evidence for coupling of the 



charge carriers to the spin one resonance at 
41 meV. The other curves in the figure are 
equally important. The gray squares repre- 
sent the Px{^) used in calculations displaced 
in energy by the gap Aq = 27 meV. It is con- 
structed completely from experiment. The 
fitting procedure involves two critical inde- 
pendent steps. First the data on T~p^{uj) in 
the normal-state is used to get a background 
spectrum of the form given in Eq. (^) which 
applies to the normal-state. This defines ujsf 
and the corresponding P. This spectrum 
is also valid at the critical temperature Tc. 
We use this to determine the last parameter, 
the anisotropy parameter g in Eqs. ( |10a|) and 
( |13a|) , in solving the linearized Eqs. ([10|) for 
g to give the required T^.. In the second step 
the normal-state result for the spectral den- 
sity is modified only in the region of the res- 
onance peak leaving it unchanged at higher 
energies. The resonance is positioned and 
its magnitude given by the data for the ex- 
perimental W{uj)/2 (solid curve). There is 
no ambiguity and the procedure is definite. 
A check on the consistency of this procedure 
is then performed in calculating the theoreti- 
cal W{uj)/2 (dashed line) from the theoretical 
optical scattering rate calculated numerically 
from the solutions of the Eliashberg equa- 
tions. We see that the theoretical data agree 
remarkably well with experiment in the re- 
gion of the resonance. More explanations of 
the differences between W{uj)/2 and Pxi}^) 
beyond the resonance region will be provided 
later on. 

After this important low temperature con- 
sistency check we can now study in more de- 
tail the temperature dependence of the spec- 
tral density Px{'^) iu the superconducting 
state of an optimally doped, twinned YBCO 
single crystal using experimental data. The 
frequency dependence of the optical scatter- 
ing rate has been studied at five tempera- 
tures by Basov et al. ||4^. We reproduce 



their experimental results in the top frame of 
Fig. p. The data are for T = 10 K solid line. 
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FIG. 6: Top frame: optical scattering rate 
T^{T,uj) in meV for optimally doped, twinned 



YBCO single crystals [42|. Bottom frame: func- 
tion W{oj) vs 00 in the region of the optical res- 
onance. 

T = 40 K dotted line, T = 60 K dashed line, 
T = 80K dash-dotted line, and T = 95K 
(in the normal-state) gray solid line. We see 
that in the normal-state t~^{(jj) vs uj shows a 
quasi-linear dependence on uj but in all other 
curves a depression below the normal curve is 
seen at small uj below roughly 75meV. The 
depression is the more pronounced the lower 
the temperature. At higher frequencies all 
curves roughly coincide. At low tempera- 
tures, {uj) as a function of u shows a small 
value up to 50meV or so, with a sharp rise 
around 75 meV characteristic of the existence 
of a sharp peak in Pxi}^)- K is clear that this 
peak increases in strength as T is lowered into 
the superconducting state. More quantitative 



information on the temperature variation of 
the optical resonance, its strength in J^x(u;), 
and its position is shown in the bottom frame 
of Fig. ^ where we show the second derivative 
W{uj) function derived directly from the data 
given in the top frame by performing the dif- 
ferentiation indicated in Eq. (^). The solid 
curve is for T = 10 K, the dotted curve for 
T = 40 K, the dashed for T = 60 K, and dash- 
dotted for T = 80 K. The height of the peak 
of the resonance clearly increases with low- 
ering of the temperature. In the curves for 
W{uj) vs UJ the position of the peak is seen 
to be reduced as T is increased towards T^. 
If it is remembered that the peak in W{(jj) 
is located at the gap value Ao(T) plus the 
position of the resonance and the tem- 
perature dependence of the gap is accounted 
for (it decreases with increasing temperature 
and rapidly goes to zero as is approached) 
then we conclude that the position of the res- 
onance is temperature independent although 
its strength decreases as T increases. This is 
also in agreement with the observation made 
by Dai et al. |^ that the energy at which 
the spin one resonance is observed in YBCO 
(41meV) is temperature independent. 

In Fig. we show results for the spectral 
density /^x(co') vs uj obtained from the data of 
Fig. 1^ at five temperatures, namely T = 95 K 
(solid gray curve), T = 80K (dash-dotted 
curve), T = 60K (dashed curve), T = 40K 
(dotted curve), and T = lOK (solid black 
curve). The procedure follows in all cases the 
procedure already described in detail for the 
T = 10 K data. Now, the resonance is posi- 
tioned and its magnitude given by the W{uj) 
data shown in the bottom frame of Fig. ^ and 
W{uj)/2 is used to modify the normal-state 
MMP background spectrum only in the re- 
gion of the resonance. A consistency check is 
then performed for all temperatures and we 
present the T = 40 K result in Fig. p. The top 
frame gives W{uj)/2 obtained directly from 
experiment in the resonance region (solid line, 
corresponds to the dotted line in the bottom 
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FIG. 7: The charge carrier-spin excitation spec- 
tral density I'^xi^) determined from optical scat- 
tering data at various temperatures. Solid gray 
curve T = 95 K, dash-dotted T = 80 K, dashed 
T = 60 K, dotted T = 40 K, and black solid 
T = 10 K. Note the growth in strength of the 
41 meV optical resonance as the temperature is 
lowered. 

frame of Fig. The model Px{^) fo^^ this 
temperature is given by the gray solid squares 
displaced in energy by the superconducting 
gap Aq = 21 meV. This spectrum agrees with 
the experimental W{u!)/2 in the resonance re- 
gion. The dashed curve, finally, represents 
the theoretical W{u)/2 and we recognize, 
again, a remarkable agreement between the- 
ory and experiment in the appropriate energy 
range. The fit to the optical scattering rate 
which in the end is the quantity that matters 
is shown in the bottom frame of Fig. ^j. The 
theoretical curve (dashed line) follows well 
the experimental data (solid line). In this 
sense we have found a spectral density which 
can reproduce the measured optical scatter- 
ing rate at T = 40 K and it does not matter 
much how we arrived at our final model for 
Px{^) which involved, ClS db step, consider- 
ation of the function W{uj) which served to 
guide our choice. 

We return now to a more detailed discus- 
sion of the optical conductivity. In Fig. || we 
show the real part of the optical conductivity 
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FIG. 8: The top frame gives the resonance 
peak obtained directly from the optical data for 
YBCO in the superconducting state (solid curve) 
at T = 40 K, the model I'^xi^) (gray squares) 
displaced by the gap, and the theoretical results 
for W{uj)/2 obtained from the model spectral 
density (dashed line). The bottom frame gives 
the optical scattering rate r^p^^uj) vs m. The ex- 
perimental results give the solid curve and our 
theoretical fit to it is the dashed curve. 



(Ti (u) for untwinned, optimally doped YBCO 
single crystals (solid line) reported by Homes 
et al. and compare with results of vari- 
ous calculations. This means that we are now 
comparing theoretical predictions with exper- 
imental YBCO data which have not been 
used to derive the Px{^) spectra. The gray 
solid curve is the BCS result for a gap of 
24meV and impurity scattering in Born ap- 
proximation corresponding to = 0.32 meV. 
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FIG. 9: Comparison of the real part of the in- 
plane conductivity cri{uj) vs uj for various models 
(T = 10 K). The grayed solid line with a peak 
before 50meV is BCS. The dash-dotted line is 
an Eliashberg calculation with an MMP spectral 
density peaked at ujsf = 20meV. The dashed 
line is the same but with our temperature and 
frequency dependent Pxi'^) (see Fig. |^ used in- 
stead of the MMP model. As described in the 
text this charge carrier-exchange boson interac- 
tion spectral density I'^xi^) has been determined 
through a consideration of in-plane optical data. 
The dotted (Born) and dash-double-dotted (uni- 
tary scattering) curves include impurities in ad- 
dition to the Pxi'^) model for inelastic scatter- 
ing. The solid line is the data of Homes et al. 



It is very clear that no agreement with experi- 
ment is possible with BCS c?-wave. One needs 
to go to an Eliashberg formulation if one is 
to even get close to the data. In some sense 
this is very positive since a good fit with a 
BCS formulation would mean that details of 
the pairing potential do not play an impor- 
tant role in the conductivity and, thus, our 
procedure would not be a good way to pin 



down some of the details of the pairing in- 
teraction. The dash-dotted curve represents 
results of Eliashberg calculations but with an 
MMP model for the Px{'^) kernel. This ig- 
nores the growth of the 41 meV optical reso- 
nance that enters when the superconducting 
state develops. While the agreement with the 
data is good at high energies beyond 100 meV 
say, it fails completely in the low energy re- 
gion. In particular, the Drude like peak at 
very low energies is much too narrow. It is im- 
portant, however, to emphasize the difference 
between Eliashberg and BCS at high ener- 
gies where BCS gives a conductivity which is 
much too small while, in comparison, Eliash- 
berg with an MMP kernel gives larger values 
reflecting the long tails extending to 400 meV 
in the MMP Px{'^) spectrum. This is taken 
as strong evidence for the existence of long 
tails in the pairing spectral density and ar- 
gues against a pure phonon mechanism which 
would be lot more confined in energy. 

The dashed curve in Fig. ^ gives our re- 
sults for the real part of the conductivity 
aiiuj) when our J^x(a;) at T = 10 K (solid 
line in Fig. |^) is used in the calculations 
rather than the MMP kernel. The existence 
of the 41 meV resonance shifts the large rise 
in the conductivity which now begins at much 
higher energies ~ 50meV than in the MMP 
case. It also leads to a maximum around 
100 meV in good agreement with experiment. 
Even better agreement can be obtained if a 
small amount of impurity scattering is in- 
cluded within the unitary or resonant scat- 
tering limit c — in Eq. (|l^). Results with 
F"*" = 0.63 meV are shown as the dash-double- 
dotted curve which displays all the impor- 
tant characteristics observed in the experi- 
mental data. The agreement is truly very 
good. A final curve including only Born im- 
purity scattering (dotted curve) shows that 
this limit cannot explain the data. It should 
be clear from this comparison that BCS the- 
ory is quite inadequate in describing the ob- 
served features of the real part of the infrared 
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FIG. 10: The imaginary part of the conduc- 
tivity a;cT2(w) vs uj for the various models de- 
scribed in Fig. ^ (T = 10 K). The sohd curve is 
the data [48|. The dash-dotted curve is the re- 
sult of an Eliashberg calculation with an MMP 
model while the other curves are based on the 
model Pxi^) which includes the 41 meV res- 
onance. These three curves are for the pure 
case (only inelastic scattering, dashed line), the 
others are with some additional elastic impurity 
scattering in Born (dotted) and unitary (dash- 
double-dotted) limit with = 0.32 meV and 
= 0.63 meV respectively. The gray solid line 
is the BCS result. 

conductivity as a function of frequency while 
Eliashberg theory can give a good fit. It 
is also clear that some impurity scattering 
in the unitary limit is needed and that the 
electron-boson exchange spectral density has 
long tails extending to at least 400 meV, and 
at T = 10 K has a large contribution from the 
41 meV resonance peak. 



In Fig. |10|we show our results for the imag- 
inary part of the optical conductivity (72 (tu). 
What is plotted is uja2{uj) vs u; at T = 10 K. 
The solid curve gives data from Homes et al. 
48| . The dash-dotted curve are Eliashberg 



results based on an MMP kernel which does 
not account for charge carrier coupling to the 
41 meV resonance. It fails to reproduce the 
data. On the other hand, when the resonance 
is included in our Px{^) "we get the dashed 



curve which agrees with the data much better 
and shows a large depression in the curve lo- 
cated around 75meV in agreement with the 
data. This is the signature in 172(0;) of the 
41 meV resonance in Pxi^)- At very low 
energies the agreement is not as good. As 
for the real part of the conductivity this re- 
gion is sensitive to impurity scattering while 
at higher energies only inelastic scattering 
is really important. The dash-double-dotted 
curve includes impurities in the unitary limit 
with F"*" = 0.63 meV as before. This pro- 
duces excellent agreement with the data even 
in the small u region which is not as well de- 
scribed in the case of Born scattering (dotted 
curve with = 0.32 meV). It is clear from 
this graph that ua2{uj) has an easily identi- 
fiably signature of the spin resonance. Also, 
the spectral density extends to high energies 
and Eliashberg theory with some contribu- 
tion from unitary scattering impurities is in 
very good agreement with the data. 

While we have seen that an unmistakable 
signature of the 41 meV optical resonance ex- 
ists in the data on both real and imaginary 
part of the infrared conductivity function 
of energy uj in the T = 10 K data, the reso- 
nance is even clearer in the second derivative 
defined by W{u) of Eq. (||). The near equal- 
ity between W{uj) and a'^F{u) established by 
Marsiglio et al. for a phonon mechanism 
was only for the normal-state, i.e.: T^p^{uj) 
entering the formula is the normal-state op- 
tical scattering rate. But in the high 
cuprates the only low temperature data avail- 
able is often in the superconducting state. 
Thus, we need to discuss in more detail what 
happens in Eq. (^ when T~p^{u) is replaced 
by its superconducting state value. As we 
have seen in Fig. |^, in the normal-state and 
at low temperatures W{uj) is almost exactly 
[] equal to the input Px{^) 
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for models based on the NAFFL. Of course, 
Px{^) is seen in W{uj) through electronic 
processes. But in the normal-state the elec- 
tronic density of states N{e) is constant in the 



important energy region and, thus, does not 
lead to additional structures in W^u) that 
are not in Pxi^)- Such additional struc- 
tures would then corrupt the signal, if the aim 
is to obtain Px{^) from W{uj). This is no 
longer the case in the superconducting state 
because of the logarithmic van Hove singu- 
larities in N{e) and these do indeed strongly 
influence the shape of W{uj) and introduce 
additional structures in W{u) corresponding 
to combinations of the positions of the singu- 
larity in N{6) and the peak in I^xi^) -^r 
as described by Abanov et al. [^] The struc- 
tures in W{uj) corresponding to these singu- 
larities contaminate the signal in the sense 
that W{ijj) in the superconducting state is no 
longer equal to the input In 
fact, only the resonance peak appears clearly 
at Ao + Er and its size in W{uj) is about 
twice the value of Pxif^) that frequency. 
In some cases the tails in W{uj) also match 
well the tails in J^x(c<j). In the end, of course, 
W{uj) serves only as a guide and it is the qual- 
ity of the final fit to the conductivity data 
that determines the quality of the derived 

Nevertheless, besides giving a measure of 
the coupling of the charge carriers to the op- 
tical resonance W{uj) can also be used to 
see the position of density of states singular- 
ities, as shown in Fig. ^ where Px{^) (gray 
squares) and W{ijj) (solid line) derived from 
our theoretical results are compared. Also 
shown by vertical arrows are the positions of 
Ao + Er, 2Ao + Er, Ao + 2Er, and 2Ao + 2Er. 
We note structures at each of these places 
and this information is valuable. Note that 
at 2Ao -|- Er the large negative oscillation 
seen in W{uj) is mainly caused by the kink 
in Px{^) (gray squares) at about 75meV. 
The density of electronic states effects clearly 
distorted the spectrum above the resonance 
peak and W{uj) stops agreeing with the in- 
put Px{^) ill this region until about 150 meV 
where agreement is recovered. In summary, 
W{uj) contains some information on singu- 
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FIG. 11: Second derivative VF(ll') compared 
with the input spectral density I'^xi^)- The 
41 meV peak in I^x(u;) (gray squares) is clearly 
seen in W{uj)/2 (solid line) as are the tails at 
higher energies. In the energy region between 75 
and 150 meV the van Hove singularities in the 
electronic density of states show up added on 
to Er and distort the correspondence between 
W{io)/2 and I'^xi^)- The dashed line shows 
the function W{io)/2 derived from experimental 
data. 

larities in N{e) as well as on the shape and 
size of Px{^) and, in the superconducting 
state, the two effects cannot be clearly sep- 
arated. Nevertheless, W{u!) remains a valu- 
able intermediate step in the construction of 
a charge carrier-exchange boson interaction 
spectral density from optical data. To close, 
the dashed curve in Fig. ^ is the direct exper- 
imental data for W{uj)/2 which is remarkably 
similar to theory when we consider that a sec- 
ond derivative is needed to get this curve. 

Next we want to concentrate on the out- 
of-plane conductivity (c-axis conductivity) of 
YBCO. Before presenting results we stress 
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FIG. 12: Comparison between in-plane (solid 
line) and out of plane (dashed line) real part of 
the c-axis conductivity aic{uj) vs co in an Eliash- 
berg model with our model carrier-boson spec- 
tral density I'^xi^) which includes the 41 meV 
spin resonance. The dotted curve is aic{uj) for a 
BCS d-wave model with the same gap value as 
in the Eliashberg work and is included for com- 
parison. 

that the boson exchange kernel I^xi^) is 
an in-plane quantity and is taken from our 
discussion of the in-plane conductivity. It 
is not fitted to any c-axis data. It is to 
be used unchanged to calculate the out-of- 
plane conductivity assuming first coherent 
hopping with t±{k) = t± cos^(0) in Eqs. (|T9|). 



The solid curve in Fig. |12| is the in-plane 
Eliashberg result which is included for com- 
parison with the dashed curve which is for 
the c-axis. In the boson assisted region, 
which would not exist in a BCS theory, both 
curves have a remarkably similar behavior. 
At very low frequencies, a region which comes 
mainly from the coherent delta function part 
of the carrier spectral density, and which is 
the only part included in BCS, we note a 
narrow Drude-like peak in the solid curve. 
This part is suppressed in the c-direction 
(dashed curve) because the contribution from 
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FIG. 13: Comparison between in-plane (dotted) 
and out of plane (solid) for the real part of the 
conductivity ai (u) vs uj. The data is from Homes 
et al. |8|. 



the nodal quasiparticles are effectively left 
out by the t±cos^{(j)) weighting term. Also, 
shown for comparison are BCS results for co- 
herent hopping (dotted curve) along the c- 
axis. These results show no resemblance to 
our Eliashberg results and also do not agree 
with experiment Fig. What determines 
the main rise in the region beyond the Drude 
part of the conductivity in aic{uj) are the bo- 
son assisted processes and this rise does not 
signal the value of the gap or twice the gap 
for that matter, but rather a combination of 
Aq and the resonance energy E^. 



In Fig. |I3] we compare the data of Homes 
et al. on the same graph for in-plane 

(dotted) and out-of-plane (solid) conductiv- 
ity cri(co'). It is clear that in the c-direction, 
the nodal quasiparticles seen in the dotted 
curve are strongly suppressed. This favors 
the t± cos^(2(/)) matrix element for the c-axis 
dynamics as we have just seen. Further, in 
the boson assisted region the two curves show 
almost perfect agreement with each other, 
which again favors the t_LCOs^(20) coupling 



as was illustrated in the theoretical curves 
of Fig. |l^. One difference is that the main 
rise, indicating the onset of the boson assisted 
incoherent (in-plane) processes, appears to 
have shifted slightly toward lower frequencies 
in the c-axis data as opposed to a shift to 
slightly higher frequencies in our theory. It 
should be remembered, however, that in the 
raw c-axis data, large structures appear in the 
conductivity due to direct phonon absorption 
and these need to be subtracted out, before 



data for the electronic background of Fig. [13 
can be obtained. In view of this, it is not clear 
to us how seriously we should take the rela- 
tively small disagreements that we have just 
described between theory and experiment. 

With the above reservation kept in mind 
we show in Fig. 0, a comparison of vari- 
ous theoretical results with experimental c- 
axis conductivity (black solid line). There are 
five additional curves. The black ones are ob- 
tained from an Eliashberg calculation based 
on the MMP model for Px{^) with impuri- 
ties = 0.32 meV included to simulate the 
fact that the samples used are not perfect, 
i.e.: are not completely pure, but this pa- 
rameter does not play a critical role in this 
discussion. Incoherent c-axis coupling is as- 
sumed with |Vi/Vo| = 1 (black dotted). It is 
clear that this curve does not agree well with 
the data and that the coupling along the c- 
axis cannot be dominated by incoherent hop- 
ping between planes. This is also in agree- 
ment with the results of a theoretical study 
by Dahm et al. 
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who also observed better 
agreement for coherent c-axis conductivity in 
the overdoped regime. On the other hand the 
fit with the black dashed line is good in com- 
parison. It uses the same MMP model but 
with coherent coupling of the form t_L(k) = 
tj_cos^(20). This fit may already be judged 
satisfactorily but it should be remembered 
that if we had used the model of Pxi}^) with 
the 41 meV peak included instead of MMP, 
the agreement would have deteriorated. This 
is troubling since one would expect that cou- 
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FIG. 14: Comparison with the data of Homes et 
al. for the c-axis conductivity (black solid 
curve). The theoretical curves were obtained in 
a BCS theory, solid gray (coherent), dotted gray 
(incoherent) and the others in Eliashberg theory 
with MMP model and impurities t"*" = 0.32 meV. 
The black dotted curve is for incoherent c-axis 
with |Vi/Vo| = 1, the dashed for coherent c-axis 
with t±{(j)) = t± cos^(20) with (p an angle in the 
two dimensional Cu02 Brillouin zone, and the 
dash-dotted is a fit to the data provided by a 
mixture of coherent and incoherent. We stress 
that this last fit is for illustrative purposes only, 
and is not unique. 

pling to the 41 meV spin resonance would 
be stronger in the c-direction data than it 
is in the in-plane data. This is because the 
c-axis emphasizes the hot spots around the 
antinodal directions which connect best to 
(vr, vr) in the magnetic susceptibility. This is 
the position in momentum space where this 
spin resonance is seen to be located in opti- 
mally doped YBCO. On the other hand, re- 
cent ARPES data p3|, ||, pi ISSl which fit 



well the MFL (marginal Fermi liquid) phe- 
nomenology show little in-plane anisotropy 
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for scattering around the Fermi surface and 
this is consistent with the findings here. 

The dash-dotted curve in Fig. |1^ illus- 
trates a fit to the data that can be achieved 
with a dominant coherent piece and subdom- 
inant incoherent contribution. It is not clear 
to us whether such a close fit is significant 
given the uncertainties in the data and the 
lack of uniqueness in the fitting procedure. 
It does, however, illustrate the fact that a 
small amount of incoherent c-axis hopping 
cannot be completely ruled out from consid- 
eration of the infrared data and that this data 
can be understood quite well within Eliash- 
berg theory. The last two curves (solid gray 
and dotted gray) are based on BCS cZ-wave 
theory and are reproduced here to illustrate 
the fact that such a theory is unable to ex- 
plain the c-axis data. The solid gray curve is 
with t_L(k) = t± cos^(20) and the dotted gray 
one for incoherent c-axis transport. Com- 
pared with our Eliashberg results the agree- 
ment with the data is poor. 



The Microwave Conductivity 

The microwave conductivity func- 
tion of temperature in pure single crystals of 
YBCO revealed the existence of a very large 
peak around 40 K [^] whose size and po- 
sition in temperature depends somewhat on 
the microwave frequency used. This peak has 
been widely interpreted as due to a rapid re- 
duction in the inelastic scattering below Tc 
and is generally referred to as the collapse of 
the low-temperature inelastic scattering rate. 
This has been taken as strong evidence that 
the mechanism involved is electronic in origin 
and, this fact translates in our formalism into 
the fact, that the charge carrier-exchange bo- 
son interaction spectral density Px{^) is re- 
duced at low frequencies due to the onset of 
superconductivity. We have already seen in 
Fig. the growth of the 41 meV resonance in 
Px{^) the temperature is lowered. At the 



same time the in-plane infrared optical data 
shows a gapping or at least a strong reduction 
of spectral weight at small u. This implies 
that for temperatures smaller than the char- 
acteristic energy associated with this reduc- 
tion, the inelastic scattering rates will become 
exponentially small and therefore the inelas- 
tic scattering time will become very large. 
This feature by itself will increase the mi- 
crowave conductivity. At the same time the 
normal fluid density is of course decreasing 
towards zero. This feature reduces the ab- 
sorption which is due only to the normal ex- 
citation. The two effects combine to give a 
maximum in the real part of the microwave 
conductivity at some intermediate tempera- 
ture. 

Another possible way to describe this col- 
lapse of the inelastic scattering time is the 
introduction of a temperature dependent in- 
elastic scattering time which can be modeled 
from spin fluctuation theory 

Recently Hosseini et al. provided new 
microwave data at five different frequencies 
between 1 and 75 GHz on ultra pure samples 
of YBa2Cu306.99 grown in BaZrOa crucibles. 
In Fig. IT5| we show results obtained from our 



Eliashberg solutions and compare with exper- 
iment The solid squares are the data 



of Hosseini et al., the open triangles are our 
numerical results in the clean limit, and the 
solid triangles include a small amount of im- 
purities characterized by F"*" = 0.003 meV 
and c = 0.2. The figure has five frames 
one for each of the five microwave frequen- 
cies considered, namely 1.14, 2.25, 13.4, 22.7, 
and 75.3 GHz. We see that even for these 
ultrapure crystals, results obtained without 
including impurities do not agree well with 
the data at the lowest microwave frequencies 
considered and at the lowest temperatures. 
For example, in the case of the uj = 1.14 
and 2.25 GHz runs the predicted peak is much 
to high. The agreement, however, improves 
as the frequency of the microwave probe is 
increased. More importantly, when a small 



a" ■ 





0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



J 1-5 



3, 

e" 1.0 



A AP*,^. 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

t t 



0.0 0.2 0.4 0.6 0.8 1.0 
t 



FIG. 15: Microwave conductivity ai{uj,t) in 
lO'^O^^m^^ vs the reduced temperature t = 
T/Tc for the five frequencies measured in ex- 
perimental work of Hosseini et al. [57] namely 
Q = 1.14, 2.25, 13.4, 22.7, and 75.3 GHz (bot- 
tom frame). Solid squares are experiment, open 
triangles clean limit and solid triangles inelastic 
scattering plus impurities characterized by a po- 
tential with r+ = 0.003 meV and c = 0.2. 



amount of impurity scattering with c = 0.2 
is included, good agreement is obtained in 
all cases. The same data plotted in a dif- 
ferent way shows better the dramatic im- 
provement in the agreement with the data 
when impurities are included. This is demon- 
strated in Fig. |16| where we show the data 
for the microwave conductivity (Ji{uj) vs uj 
at three different temperatures. The data 
are represented by solid squares, up-triangles 
and down-triangles for T = lOK, 15K and 
20 K respectively. The open symbols give the 
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FIG. 16: The microwave conductivity cri(a;,T) 
as a function of u for three different temper- 
atures. The data is the same as shown in 
Fig. ||. The open symbols are theory for the 
pure limit, the solid gray symbols theory with 
some impurity scattering additionally included, 
and the solid black symbols are experiments. 
The squares are for T = lOK, the up-triangles 
for T = 15 K, and the down-triangles for T = 
20 K. 



results of our Eliashberg calculations in the 
pure case and the solid gray symbols include 
impurities. The gray lines through the points 
guide for the eye. The agreement with 
the data in this last case is within experimen- 
tal error and is acceptable. It is clear, that a 
small amount of elastic scattering needs to be 
included in the calculations to achieve good 
agreement. 

We now turn to the c-axis. No new pa- 
rameters relevant to the in-plane dynamics 
need to be introduced in order to understand 
the c-axis data. It is necessary, however, to 
have some model for the c-axis charge trans- 
fer. Coherent or incoherent hopping will lead 
to quite different conclusions as will the as- 
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FIG. 17: The c-axis microwave conductivity 
o"ic(T) at w = 22 GHz as a function of tem- 
perature T. The open up-triangles were ob- 
tained from the empiricahy determined I'^xi^) 
shown in Fig. |^ The sohd squares are for an 
MMP form Eq. (jl^) with low frequency cutoff 
apphed. This cutoff is the same as seen in Fig. ^. 
The open down-triangles employ the same MMP 
model with ujsf = 20meV and without the low 
frequency cutoff. In this case there is no peak 
in aic{T). All curves are for coherent tunneling 
with t^(k) = t±cos'^{2e). 



sumption, in the coherent hopping case, of a 
constant or a momentum dependent hopping 
probability. For the constant case the con- 
ductivity will mirror its in-plane value but 
its magnitude will of course be greatly re- 
duced. For a ^-dependent matrix element, on 
the other hand, the nodal quasiparticles are 
eliminated from participation in the c-axis re- 
sponse and we can expect a behavior different 
from the in-plane results. 



In Fig. we show numerical results for 
the temperature dependence of the real part 
of the c-axis microwave conductivity o"ic(T) 
at uj = 22 GHz function of tempera- 



ture T in arbitrary units. The results are 
for coherent hopping with a ^-dependent ma- 
trix element. For the in-plane case the same 
form applies except that the vertex t^(k) 
would be replaced by a Fermi velocity fk. 
In as much as both these vertices are taken 
to be independent of k they can be pulled 
out of the integral over k in Eqs. (|r^) and 
in-plane and out-of-plane conductivities dif- 
fer only by a numerical constant which sets 
the over all scale in each case. Since no peak 
is observed in the temperature dependence 
of the c-axis conductivity this case does 
not agree with experiment and will not be 
treated further here. Only results for t_L(k) = 
t±cos'^{29) are considered in Fig. |1^. The 
open up-triangles are the results obtained 
from the charge carrier-exchange boson inter- 
action spectral density Px{^) obtained em- 
pirically from the in-plane infrared conductiv- 
ity (Fig. 1^). This is the only material param- 
eter which characterizes YBCO in the Eliash- 
berg equations (|T3|). Solutions of these equa- 
tions determine the in-plane Green's function 
(p^) and hence the c-axis conductivity Eqs. 
(pl|). Arbitrary units are used, so that the 
absolute value of t± is not required. We see 
a broad peak in aic{T) vs T which is cen- 
tered around T = 60 K rather than around 
T = 40K for the in-plane case of Fig. 15. 
The c-axis peak is also smaller. These dif- 
ferences are entirely due to the extra factor 
of cos^(2^) in the c-axis conductivity which 
eliminates the nodal direction. This has a 
profound effect on the resulting temperature 
dependence of cric(T) but, as we can see, 
does not entirely eliminate the peak in aic{T). 
There is considerable disagreement with ex- 
periment which is shown as the solid squares 
in Fig. |1^. Our theoretical results are ro- 
bust in the sense that the peak is due to 
the greatly reduced spectral weight in J^x(co') 
of Fig. 1^ at small u when superconductiv- 
ity sets in and this is fixed from consider- 
ation of the in-plane conductivity. The ef- 
fect of this spectral weight reduction is fur- 
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FIG. 18: The c-axis microwave conductivity 
(Tic(T) at uj = 22 GHz as a function of tempera- 
ture T. The sohd squares are the experimental 
results by Hosseini et al. (5^ shown for com- 
parison. The others are theory based on various 
models for the charge carrier-exchange boson in- 
teraction spectral density. All calculations are 
for the incoherent case based on Eqs. (^) with 
^1 = ^ in ( p2|) and Vq adjusted to match the ex- 
perimental value at T = 90 K. Open up-triangles 
are based on our empirically determined charge 
carrier-exchange boson interaction spectral den- 
sity I'^xi^) of Fig. without additional impurity 
scattering (pure limit). The open down-triangles 
are also in the pure limit but the MMP model 
( p^ ) is used without cutoff. The stars are the 
same as the open down-triangles but now impu- 
rity scattering is included in the unitary limit 
with r+ = 0.5 meV. 



ther illustrated in Fig. by the solid squares 
which employ instead of our empirical value 
for /^x(^) the simpler MMP form of Eq. ( p4| ) 
with the same low frequency cutoff as indi- 
cated in Fig. 1^ being applied. The cutoff is 
of course temperature dependent and goes to 
zero at Tc. The peak in cric(T) vs T remains 



and is close to the results obtained when an 
optical resonance is included in addition to 
a low frequency cutoff. For comparison, the 
down-triangles were obtained when no low 
frequency cutoff was applied. We now see 
that the peak in aidT) vs T is completely 
eliminated. This demonstrates that the peak 
is due to the collapse of the inelastic scatter- 
ing rate embodied in the low frequency gap- 
ping of the charge carrier-boson spectral den- 
sity. In summary, even when a momentum 
dependent coherent hopping matrix element 
of the form t_L(k) = t± cos^(26') is considered, 
gaping of J^x(^) at small u leads directly 
to a peak in the c-axis microwave conduc- 
tivity. However, the spectral density Px{^) 
which enters the Eliashberg equations ([T0|) on 



the imaginary axis and (|T3D on the real axis 
could depend on position on the Fermi sur- 
face. This complication was not considered 
here but it is important to point out that 
coherent c-axis tunneling could lead to rea- 
sonable agreement with the measured tem- 
perature variation of the microwave conduc- 
tivity, if the spectral density Px{^) is dif- 
ferent along the antinodal direction and, in 
particular, has no gapping at low frequencies. 
This difference might not have been picked 
up in our analysis of the in-plane conductiv- 
ity which is characteristic of an average over 
all points on the Fermi surface and not just 
of the antinodal direction. 



Fig. |Tg shows results for the temperature 
variation of the c-axis microwave conductiv- 
ity o"ic(^) at u; = 22 GHz in the incoher- 
ent coupling case, Eqs. (^ll) . As we have 
indicated in Sec. 2 this formula involves a 
double integral over momentum which sepa- 
rately weights the two Green's functions. For 
simplicity, we show results only for the case 
Vi = Vq in the impurity model potential of 
Eq. (p^ . Other values have been consid- 



ered but this does not change qualitatively 
any of the conclusions we will make. The 
open up-triangles give results when the em- 
pirical Px{^) of Fig- is used. We see that 
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in this case the theory predicts no peak in 
(Jic(T) vs T in good agreement with the ex- 
perimental results of Hosseini et al. (solid 
squares) . The low frequency cutoff built into 
our Pxiy^) (open up-triangles) has little ef- 
fect on the resulting aic{T). This is veri- 
fied directly when we compare with the open 
down-triangles which were obtained with the 
MMP form (^ without cutoff. These re- 
sults differ very little from the previous ones 
and show that the application of a low fre- 
quency cutoff does not play a critical role for 
the incoherent case. This is in sharp con- 
trast to the coherent case in which the low 
frequency cutoff leads directly to a peak in 
aiciT). The curves are also robust to the in- 
troduction of some elastic impurity scatter- 
ing as is demonstrated with the final set of 
results in Fig. |1^, denoted by stars, which is 
based on an MMP model with elastic impu- 
rity scattering included in the unitary limit 
with r+ = 0.5 meV in Eq. (0). We see that 
the inclusion of impurities does not appre- 
ciably change our results. The calculation 
clearly shows that the observed data can be 
understood naturally in an incoherent c-axis 
transport model and that the results are ro- 
bust to changes in cutoff at low uj and to the 
addition of impurities. 

This contrasts with the case of the c-axis 
infrared data which we described previously 
and found to support coherent rather than 
incoherent c-axis charge transfer. 



Other Superconducting State Properties 

The temperature dependence of the area 
under the spin resonance seen at (vr, vr) 
by spin polarized neutron scattering |^ 
has been measured and its temperature de- 
pendence denoted by {m^^^iT)) / {ml^^{T = 
10 K)) is reproduced in the top frame of 
Fig. |19| as the solid circles. Also shown on the 
same plot are our results for the area under 
the optical resonance in our spectral density 
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FIG. 19: Top frame: spectral weight under 
the optical resonance as a function of tempera- 
ture (solid line) obtained from the optical data of 
Fig. ^ (top frame, five temperatures only). The 
solid circles are the data of Dai et al. for 
the normalized area under the spin resonance 
obtained by neutron scattering. The dashed 
curve gives our calculated Ao(T)/Ao(0). Bot- 
tom frame: the normalized London penetration 
depth squared (A(0)/A(T))^ vs reduced temper- 
ature t = T/Tc (solid line) compared with the 
experimental results of Bonn et al. Q. The 
dotted curve gives thermodynamic critical field 
H,{T)/H,{0) vs t. 



Px{^) various temperatures (see Fig. |^). 
We denote this by A(T) and plot as the solid 
line the ratio A{T)/A{T = 10 K) which fol- 
lows the same temperature variation as the 
neutron result. This temperature variation is 
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FIG. 20: The electronic part of the thermal 
conductivity as a function of the reduced tem- 
perature t = T/Tc. The solid squares are the 
experimental results of Matsukawa et al. |60| 
and the solid triangles our theoretical result with 
some impurity scattering included. 



also close to that of the gap edge shown as 
the dashed curve. This last curve was found 
byproduct of our Eliashberg calculations 
based on the numerical solutions of Eqs. (|I^). 
The gap amplitude at temperature T is given 
by 3fJe{A(w = Ao;T)} = Ao(T). The same 
solutions give the temperature dependence of 
the penetration depth and of the thermody- 
namic critical field which we present in the 
bottom frame of Fig. The solid curve 
gives resuhs for [A(0)/A(T)]2 vs t = T/T^ (the 
reduced temperature) which are close to the 
experimental results of Bonn et al. Q given 
as solid squares. 

Further results are shown in Fig. |20| for the 
electronic part of the thermal conductivity 
ClS db function of temperature. It shows, as 
does the in-plane microwave conductivity, a 
large peak around 40 K. The solid triangles 
are the results of our calculations while the 
solid squares are the experimental results of 
Matsukawa et al. [Q. The agreement is re- 
markably good. We stress that no adjustable 
parameters enter our calculations except for a 
choice for the impurity parameter t"*" which is 
also restricted in this particular case because 



TABLE I: Some superconducting properties of 
the twinned YBCO sample: AF{0) is the con- 
densation energy at T = in meV/Cu-atom, 
Us/n is the superfluid to total carrier density ra- 
tio, Qp is the plasma frequency in eV. 

Theory Experiment Ref. 



AF(0) 

ris/n 

2Ao/kBT, 



0.287 
0.33 
2.36 
5.1 



0.25 
0.25 
2.648 
5.0 
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t~^ has already been determined by a previous 



fit to the microwave data [|T^. The kernel in 
the Eliashberg equations is completely deter- 
mined from optical data. 

Other important successes of our Eliash- 
berg calculations are summarized in Table |. 
We begin with a discussion of the plasma 
frequency Qp. Referring to Fig. ^ we point 
out the arrow which shows the frequency at 
which we made our calculated conductivity 
agree exactly with experiment. This sets the 
plasma frequency which is also the total spec- 
tral weight under the real part of the conduc- 
tivity. The optical spectrum sum rule is 
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A value of fip = 2.36 eV is found which agrees 
well with the experimental value 2.648 eV (see 
Tab. 1). A further comparison of our model 
with the infrared data is provided by the 
analysis of the fraction of the total normal- 
state spectral weight which condenses into 
the superfluid: Us/n. Indeed, strong electron- 
boson coupling reduces the spectral weight 
of the quasiparticle component of the elec- 
tronic spectral function A{k,uj) compared to 
its non-interacting value by a factor of Z lead- 
ing at the same time to the appearance of an 
incoherent component. It is the latter com- 
ponent which is responsible for the Holstein 



band in the optical conductivity whereas the 
coherent quasiparticle part gives rise to the 
Drude term at T > Tc and to the superfluid 
density at T = in the spectra of a\ [uS) . 
The values of n^/n and hence {Z — 1) yield 
an estimate of the strength of renormaliza- 
tion effects in the interacting system. Tanner 
et al. 1^ obtained n^/n ~ 0.25 in crystals 
of YBCO. This compares well with the value 
~ 0.33 which corresponds to Z ~ 3 (at low 
temperatures) generated in our analysis. 
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We have also calculated the condensation 
energy as a function of temperature. Its 
value at T = follows from the normal- 
state electronic density of states which we 
take from band structure theory equal to 
2.0 states/eV/Cu-atom (double spin) around 
the middle of the calculated range of val- 
ues |7^. This gives a condensation en- 
ergy AF(0) = 0.287 meV/Cu - atom which 
agrees well with the value quoted by Norman 
et al. from the work by Loram et al. . 
(See Tab. |.) This is equivalent to a thermo- 
dynamic critical field fioHc{0) = 1.41 T with 
Hc{T) defined through AF(T) = Hl{T)/M. 
The normalized value Hc{T)/Hc{0) is shown 
as the dotted line in the bottom frame of 
Fig. and is seen to follow reasonably, but 
not exactly, the T dependence of the normal- 
ized penetration depth. One further quan- 
tity is the ratio of the gap amplitude to the 
critical temperature which in BCS theory is 
2Aq/ ksTc = 4.2 for (i-wave. In Eliashberg 
theory the gap depends on frequency. In this 
case an unambiguous definition of what is 
meant by Aq is to use the position in energy of 
the peak in the quasiparticle density of states 
which is how the gap Aq is usually defined 
experimentally for a (i-wave superconductor. 
We get a theoretical value of 2Aq/ ksTc ~ 5.1 
in good agreement with experiment, as shown 
in Tab. |. 
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FIG. 21: The temperature dependent a6-plane 
optical scattering rate for an optimally doped 
Bi2212 single crystal with E\\a in the normal 
state. The solid lines represent experimental 
data hy Tn et al. \72\. The dotted lines give 
the theoretical result calculated using Eq. @ 
and a I'^xi^) which is just an MMP spectrum 
with uJsF = 82meV. At T = 295 K theory repro- 
duces experiment almost ideally; this agreement 
deteriorates at lower temperatures. Finally, the 
dashed lines present theoretical results found for 
the same MMP spectrum as before but now with 
coupling to the 43 meV optical resonance added. 

The Compound Bi2Sr2CaCu208+5 

Optical data published by Puchkov et al. 
|46[| for optimally doped samples of the com- 
pound Bi2Sr2CaCu208+5 (Bi2212) have first 
been analyzed by Schachinger and Carbotte 
|45[| . They reported that the normal-state 
optical scattering rate (T = 300 K) can be 
fitted perfectly by an MMP spectrum with 
uJsF = 100 meV and with an high energy cut- 
off at 400 meV. The inversion of the supercon- 
ducting state optical scattering rate revealed 



the coupling of the charge carriers to a res- 
onance found at an energy of 43meV. This 
corresponds to the magnetic resonant mode 



found by Fong et al. [|7T| using inelastic neu- 
tron scattering. This mode appears below 
Tc and its intensity increases with decreasing 
temperature. 

Tu et al. recently studied the a6-plane 
charge dynamics in optimally doped Bi2212 
single crystals (T^, = 91 K) using an exper- 
imental technique with much improved sig- 
nal to noise ratio. They developed an ex- 
perimentally unambiguous method which ex- 
amines the maxima and minima of W{uj), 
Eq. (H). The authors argued that a compar- 
ison of their spectral data with data found 
for YBCO suggests that a pseudogap exists 
in Bi2212 above T,, at least at T = 100 K. 



Fig. ^ presents their data (solid lines) to- 
gether with a theoretical analysis. For Pxi'^) 
an MMP spectrum is used and at T = 295 K 
our best fit is found for ujsf = 82meV to- 
gether with an high energy cutoff 400 meV us- 
ing theoretical results calculated from Eq. @ 
(dotted line). If we calculate the optical 
scattering rate for the temperatures T = 
200 and 100 K using the same Px{^) spec- 
trum it becomes obvious that the agreement 
with experiment deteriorates with decreas- 
ing temperature (dotted lines). A compar- 
ison of the anomaly in the optical data at 
T = 100 K around 50meV with the optical 
data of YBCO (Fig. |) reveals quite simi- 
lar behavior which suggests that in Bi2212 
a coupling to an optical resonance can actu- 
ally be seen in the normal state. Indeed, the 
quality of the data is good enough to allow 
us to derive W{u) by inversion. The result 
shows a pronounced peak at 43 meV which 
can be used to modify the Px{^) spectrum. 
This is shown in Fig. Using this mod- 
ified spectrum Px{^) to calculate the opti- 
cal scattering rate for T = 100 K results in 
excellent agreement between experiment and 
theory (dashed line, Fig. ^ labeled 100 K). 
The T = 200 K data show a similar, but 
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FIG. 22: The charge carrier-spin excitation spec- 
tral density /^x(a;) determined from normal 
state optical scattering data shown in Fig. 21 
for optimally doped Bi2212 single crystals. The 
solid curve is for T = 295 K, the dotted curve for 
200 K, and the dashed one for 100 K. Note the 
growth in strength of the 43meV optical reso- 
nance as the temperature is lowered. 

less pronounced, anomaly. It is not possible 
to derive a W{uj) directly from the data so we 
simply use the Px{^) found for 100 K and re- 
duce the size of the 43meV peak until best 
agreement between experiment and theory is 



reached (dashed line in Fig. 21 labeled 200 K). 
This results in an Px{^) presented in Fig. |22| 
which still contains a pronounced contribu- 
tion from the coupling of the charge carriers 
to the optical resonance (dotted line). From 
this we can conclude that the optical reso- 
nance, and probably connected with it, the 
magnetic resonant mode exists at least up to 
200 K. Nevertheless, our result for 295 K in- 
dicates that at this temperature the optical 
resonance no longer exists. 

We will now concentrate on the supercon- 
ducting state and study the temperature de- 
pendence of Px{^) below Tc. The top frame 
of Fig. p3| presents the infrared scattering rate 
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FIG. 23: Top frame: optical scattering rate 
T~p{T,uj) in meV for optimally doped, un- 
twinned Bi2212 single crystals ||7^. Bottom 
frame: function W{ij)/2 vs uj in the region of 
the optical resonance. 



as measured by Tu et al. [|72| for three tem- 
peratures, namely 6K (gray solid line), 50 K 
(dashed line), and 80 K (dotted line). In com- 
parison with similar results for YBCO (top 
frame of Fig. |^) we recognize that even at 
80 K Bi2212 shows a very strong suppression 
of T~p{uj) at energies below 50meV which 
is an indication of stronger coupling of the 
charge carriers to the optical resonance. The 
bottom frame of this figure shows the func- 
tion W{uj)/2 derived from experiment with 
the high energy negative parts suppressed be- 
cause we want to concentrate on the optical 
resonance. It increases as T is lowered and 
shows only little further variation below 50 K. 
In W{ijj)/2 the resonance peak is positioned 



at the resonance energy plus the gap value 
Ao(T) and with the temperature dependence 
of the gap accounted for, we can conclude 
that the position of the resonance is tempera- 
ture independent and stays at Er = 43 meV, 
the energy at which the magnetic resonant 
mode is found by inelastic neutron scatter- 
ing ||7T|. The coupling of the charge carriers 
to a boson at 43 meV has also been observed 
in photoemission 



[3[ and tunneling work 
on Bi2212 and the 43 meV magnetic resonant 
mode seems to be the obvious candidate for 
the origin of this boson. 

Fig. 123 demonstrates the agreement which 



can be achieved between theory and experi- 
ment. The top frame shows the infrared scat- 
tering rate T^^(a;) vs a; for T = 6K. The 
solid line gives the experimental data while 
the dashed and dotted curves represent theo- 
retical results for a clean limit system and for 
a system with impurity scattering in the uni- 
tary limit (F"*" = 0.5 meV) respectively. The 
differences in the scattering rate are marginal 
for these two model systems, nevertheless, 
they become important when the real part 
of the optical conductivity cri(ci;) is investi- 
gated. The bottom frame of Fig. |23 shows 



the results. The solid line is experiment, the 
dotted line is theory for the clean limit sys- 
tem. It reproduces nicely the maximum in 
(Ti{u!) around 120 meV and the high energy 
tail. At energies below 75 meV the clean limit 
results deviates strongly from experiment to- 
wards nearly zero values and show a very pro- 
nounced, narrow peak around uj = 0. Re- 
sults for the system with impurities treated 
in the unitary limit (F+ = 0.5 meV, dashed 
line) display all important features observed 
in the experimental data. We see, as in the 
case of YBCO, that impurities affect only the 
low energy region {u < 60meV), the region 
60 < cij < 100 meV is dominated by the cou- 
pling to the optical resonance modeled in the 
Px{^) while the energy region u; > 120 meV 
is determined by the normal-state MMP part 
of /^x(co') as has already been described for 
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the YBCO compound. 

The bottom frame of Fig. ^ contains 
an arrow which points out the frequency at 
which we made our calculated cri(a;) to agree 
exactly with experiment. (This was only 
done for the clean limit calculation, the same 
scaling was used for the system with im- 
purities.) This sets the plasma frequency 
Qp = 2.3 eV which is to be compared with 
the Qp = 1.98 eV used by Tu et al. 0. Fi- 
nally, we found for the superconducting gap 
at T = 6 K a value of 25 meV. This is cer- 
tainly smaller than the value of 34meV re- 
ported by Riibhausen et al. from Raman 
spectroscopy on Bi2212 single crystals with a 
T, of 95 K. 



Application to Other Cuprates 

In contrast to the systems studied so far 
Tl2Ba2Cu06+<5 (T12201) is a monolayer com- 
pound while YBCO, Bi2212 and YBa2Cu408 
(Y124) are bilayer compounds. Moreover, 
T12201 is the only system with tetragonal 
symmetry, all the other compounds are of 
orthorhombic symmetry. Its Tc ~ 90 K and 
this is similar to the of YBCO and Bi2212. 
Y124, on the other hand has a slightly lower 
Tc of 82 K and shows properties which resem- 
ble a moderately underdoped YBCO com- 
pound. 
In Fig. 
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normal-state r^^iuj 



we show our result for the 
related to the con- 
ductivity by Eq. (0) for T12201 with = 
90 K at temperature T = 300 K. The solid 

~ The 



H 



curve is the data of Puchkov et al. 
dotted curve is our best fit for o^^i? = 30 meV 
with P adjusted to get the correct absolute 
value of the scattering rate at T = 300 K and 
uj = 200 meV. We see that this value of usf 
does not give a satisfactory fit to the data. 
The dash-dotted curve, however, fits the data 
well and corresponds to usf = 100 meV. This 
fit provides us with a model Px{^) valid for 
the normal state of T12201. This Px{^) is 



then used to calculate the anisotropy parame- 
ter g from the solution of the linearized imag- 
inary axis Eliashberg equations (|1^) for the 
critical temperature Tc = 90 K. As a result of 
this procedure all necessary parameters are 
fixed and we can now proceed to study the 
superconducting state. 

Results are shown in Fig. The solid 
line in the top frame shows the optical res- 
onance obtained from inversion of the ex- 
perimental superconducting optical scatter- 
ing rate, presented in the bottom frame of 
this figure (solid line). The gray squares are 
the Px{^) used in the calculations displaced 
in energy by the gap Aq = 26 meV. It is con- 
structed completely from experiment and we 
followed the procedure already described in 
detail for the YBCO compound. The dashed 
curve, finally, is the result of an inversion of 
theoretical data, shown in the bottom frame 
of this figure (dashed line) and we see that it 
agrees reasonably well with experiment (solid 
line). These results allowed Schachinger and 



Carbotte to predict for T12201 a spin res- 
onance at 43meV. They also predicted that 
the resonance should be less pronounced and 
broader in T12201 than in YBCO or Bi2212. 
Recently He et al. [|7^] succeeded in prepar- 
ing a T12201 sample big enough for inelas- 
tic magnetic neutron scattering. This sam- 
ple consists of about 300 coaligned optimally 
doped T12201 single crystals. This experi- 
ment confirmed the existence of a magnetic 
resonant mode in T12201 below Tc which is lo- 
cated at about 47 meV and which appears to 
be narrower than the resonances observed in 
YBCO or Bi2212. This is in slight disagree- 
ment with the results of Schachinger and Car- 
botte ||4^ and this disagreement could proba- 
bly be explained by the poorer quality of the 

m 



for the 



samples used by Puchkov et al. 
optical measurements many years ago. Other 
optical data are not available. Nevertheless, 
the basic agreement between the observation 
of Schachinger and Carbotte ||^ that the 



charge carriers in T12201 couple to an op- 
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tical resonance and the subsequent observa- 
tion of a magnetic resonant mode at about 
the same energy by He et al. |]76[ using in- 
elastic neutron scattering is quite important. 
It proves that magnetic resonant modes are 
not restricted to bilayer compounds and that 
we seem to be confronted with a unified phe- 
nomenological picture. 

The optical resonance peak is not observed 



in all systems as is illustrated in Fig. for an 
overdoped sample of T12201 with Tc = 23 K. 
In this case a fit to the T = 300 K normal- 
state data (solid gray curve) with an MMP 
model gives usf = 100 meV (gray dashed 
curve). The same spectrum also produces a 
good fit (black dashed line) to the data at 
T = 10 K (black solid line) in the super- 
conducting state. There is no need to in- 
troduce a spin resonance. Indeed the black 
solid curve for the measured optical scatter- 
ing rate T~p^{u) is smooth and increases grad- 
ually as u increases with no clear sharp rise 
at any definite frequency in sharp contrast 



with Fig. |2g. We conclude from this analysis 
that the resonance observed in some cuprates 
with high values of Tc at optimum doping is 
not present in all cases and in particular there 
is no evidence for such a resonance in over- 
doped T12201 with Tc = 23 K. In this case 
a standard MMP spectrum of the form ( p4| ) 
gives an adequate representation of the su- 
perconducting state optical scattering rate as 
a function of uj with the same spectral density 
as was determined by the data at T = 300 K. 
This is in contrast to the other cases studied 
above for which the onset of superconductiv- 
ity appears to produce essential modifications 
of the underlying spectral density I^xi^)- 

We extend our analysis to the material 
Y124 (Tc = 82 K) where we predict from 
Fig. ^ (top frame, solid curve) a resonance at 
38meV. (This is below the energy of 41meV 
for the resonance in YBCO which is not 
surprising as it is a well established prop- 
erty of bilayer high- Tc cuprates that the en- 
ergy of the magnetic resonant mode tracks 



Tc in underdoped systems and the op- 
tical resonance seems to be closely related 
to this magnetic resonant mode.) The top 
frame of this figure demonstrates the agree- 
ment with W{uj)/2 and Px{^) which was 
shifted by the theoretical gap Aq = 24meV 
which is a prediction of our calculations as, 
to our knowledge, no experimental data ex- 
ist for this material. The bottom frame of 
Fig. ^ presents our comparison between ex- 
perimental and theoretical optical scattering 
rates. The normal-state scattering rate (gray 
lines) at T = 300 K gives evidence for the ex- 
istence of a high energy background as the 
experimental data (gray solid line) are best 
fit by an MMP spin-fluctuation spectrum as 
described by Eq. (|2^) with usf = 80meV 
and a high energy cutoff of 400 meV (gray 
dashed line). The black lines compare the 
theoretical results (dashed line) with exper- 
iment [0 (solid line) in the superconduct- 
ing state at T = 10 K. The signature of the 
optical resonance, the sharp rise in 



UJ 



starting at around 50 meV is correctly repro- 
duced by theory. For u > 120 meV the exper- 
imental scattering rate shows only a weak en- 
ergy dependence and the theoretical predic- 
tion starts to deviate from experiment. This 
is in contrast to our results found for all other 
compounds and could be related to the fact 
that the Y124 compound shows features of an 
underdoped system. 



SUMMARY 

An extended Eliashberg theory can be ap- 
plied to describe the superconducting prop- 
erties of hole doped high- Tc cuprates. The 
extension goes in two directions: first, it is es- 
sential to allow the pairing potential to have 
dx2^y2 symmetry, and, second, the charge 
carrier-exchange boson interaction leading to 
pairing has to be modeled using a phe- 
nomenological approach because the micro- 
scopic origin of the attractive interaction be- 
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tween the charge carriers is still unknown. 

An anomalous steep rise in the super- 
conducting state optical scattering rate ob- 
served in optimally doped YBCO in the en- 
ergy range 50 < a; < 90 meV was attributed 
to the coupling of the charge carriers to an 
optical resonance located at about 41 meV. 
This optical resonance has its counterpart in 
a magnetic resonant mode which can be ob- 
served in YBCO at the same energy by inelas- 
tic neutron scattering. This resonance is not 
observed above Tc in the normal-state and the 
normal-state infrared scattering rate too does 
not develop any anomaly. Further experi- 
mental data on the temperature dependence 
of the infrared scattering rate and of the mag- 
netic resonance proved further agreement as 
the area under the optical resonance is seen 
to have the same temperature dependence as 
the 41 meV magnetic resonant mode. All this 
resulted in a definite procedure which allows 
a phenomenological charge carrier-exchange 
boson interaction spectral density Pxi'-^) to 
be derived which reflects the coupling of the 
charge carriers to the optical resonance and 
also describes properly the almost linear fre- 
quency dependence of the normal-state in- 
frared scattering rate. Using this phenomeno- 
logical J^x(a') as the kernel of an extended 
Eliashberg theory allows us not only to repro- 
duce the experimental infrared optical data, 
it also allows to reproduce properly the tem- 
perature dependence of the microwave con- 
ductivity, of the London penetration depth, 
and of numerous other superconducting prop- 
erties. 

This success justifies the extension of this 
analysis to other compounds, like Bi2212, 
T12201, and Y124 for which less extensive ex- 
perimental data are available. The latest high 
quahty optical data on Bi2212 proved that, 
also in this case, the coupling of the charge 
carriers to an optical resonance at 43 meV can 
be associated with an anomalous steep rise in 
the superconducting state infrared scattering 
rate. In contrast to YBCO the anomaly in 



the infrared optical scattering rate can also 
be observed in the normal state. This op- 
tical resonance has, in the superconducting 
state, its counterpart in a magnetic resonant 
mode observed by inelastic neutron scatter- 
ing at 43meV. This mode has, so far, not 
been observed in the normal state of opti- 
mally doped Bi2212. The method to derive 
a phenomenological Pxi^) from optical data 
developed for YBCO can also be applied to 
the compound Bi2212 and leads, again, to 
a temperature dependent kernel /^x(a;) and 
the extended Eliashberg theory allows an ex- 
cellent reproduction of some superconducting 
state properties. 

A similar anomaly can be observed in 
the superconducting state low-temperature 
infrared scattering rate of T12201 and Y124. 
It has, consequently, been interpreted as the 
coupling of charge carriers to an optical reso- 
nance, not present above Tc. Recently, in the 
monolayer compound T12201 a magnetic res- 
onant mode has been observed in the super- 
conducting state but at a slightly higher en- 
ergy than predicted from optical data. Never- 
theless, this fact is quite important because 
it established that the existence of a mag- 
netic resonant mode is not restricted to bi- 
layer compounds. The existence of a mag- 
netic resonant mode is still to be proved in 
Y124, in which an optical resonance seems 
to exist at an energy of 38meV. For both 
compounds, an extended Eliashberg theory 
together with a phenomenologically derived 
kernel /^x(ci;) resulted in a good reproduction 
of the optical data. Not enough information 
about a possible temperature dependence of 
the /^x(cj) in these compounds is available to 
extend the theoretical analysis to other super- 
conducting state properties. 

All this established a unified phenomeno- 
logical picture for hole doped high- Tc 
cuprates which interpretes anomalies in the 
charge carrier dynamics observed in opti- 
mally and overdoped samples as a signature 
of spin degrees of freedom in these com- 
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pounds. 
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FIG. 24: The top frame gives the optical scat- 
tering rate r^^(ti;) vs lo for optimally doped 
Bi2212 single crystals at a temperature of 6K. 
The experimental results give the solid curve and 
our clean limit theoretical fit to it is the dot- 
ted curve. The dashed curve presents theoreti- 
cal results for a system with impurity scattering 
in the unitary limit described by the parameter 
F"*" = 0.5 meV and c = 0. The lower frame gives 
a comparison of the real part of the in-plane op- 
tical conductivity ai{u;) vs lo for the two models 
already presented in the top frame. The solid 
curve is the experimental data, the dotted curve 
the clean limit theoretical result, and the dashed 
curve the theoretical result for the system with 
impurities. 
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FIG. 25: The normal-state optical scattering 
rate t-^{oj) vs uj for T12201 with a = 90K 
obtained from the work of Puchkov et al. |46| 
(solid curve) . The dash-dotted curve from theory 
based on Eq. with an MMP model spectral 
density using a spin fluctuation frequency losf = 
100 meV gives good agreement while the choice 
of 30meV (dotted curve) does not. 
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FIG. 26: The top frame gives our model for 
the spin-fluctuation spectral density (displaced 
by the theoretical gap Aq = 26meV for T12201 
in the superconducting state at T = 10 K (gray 
solid squares). The dashed line is W{u) obtained 
from the calculated conductivity and the black 
solid line is the coupling to the resonance found 
directly from experiment. (The high frequency 
part has been omitted.) It was used to in con- 
structing the model I^xi^)- The bottom frame 
shows the optical scattering rate at T = 10 K 
(solid line) and the theoretical fit to experiment 
found from Eliashberg theory. 



38 




CO (meV) 

FIG. 27: The optical scattering rates in an over- 
doped sample of T12201 with a Tc = 23 K. The 
solid lines represent experimental data and the 
dashed lines fits. The gray curves apply in the 
normal-state at T = 300 K and the black curves 
in the superconducting state at T = 10 K. No 
optical resonance peak is found in this case. 
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FIG. 28: The same as Fig. 26 but for the ma- 
terial Y124. The spin-fluctuation spectral den- 
sity Px{^) was displaced by the theoretical gap 
Aq = 24meV in the top frame. In addition, the 
grayed lines in the bottom frame of this figure 
show the comparison between experimental and 
theoretical normal-state data at T = 300 K. Due 
to this comparison ojsF = ^0 meV for the MMP 
spectrum of Eq. (24). 



